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We use the two-dimensional Kolmogorov-Smirnov (KS) test to study the correlation between the 
60 cosmic ray events above 4 x 10^^ eV from the AGASA experiment and the positions of infrared 
luminous galaxies from the IRAS PSCz catalog. These galaxies are expected to be hosts to gamma 
ray bursts (GRB) and magnetars, both of which are associated with core collapse supernovae and 
have been proposed as possible acceleration sites for ultra high energy cosmic rays. We find con- 
sistency between the models and the AGASA events to have been drawn from the same underlying 
distribution of positions on the sky with KS probabilities > 50%. Application of the same test to 
the 11 highest AGASA events above 10^" eV, however, yields a KS probability of < 0.5%, rejecting 
the models at 99.5% significance level. Taken at face value, these highest energy results suggest that 
the existing cosmic ray events above 10'^'^ eV do not owe their origin to long burst GRBs, rapidly 
rotating magnetars, or any other events associated with core collapse supernovae. The larger data 
set expected from the AUGER experiment will test whether this conclusion is real or is a statistical 
fluke that we estimate to be at the 2a level. 



I. INTRODUCTION 

The source of cosmic rays at energies above the 
Greisen-Zatsepin-Kuzmin (GZK) feature predicted to 
appear at energies above ~ 4 x 10^^ eV, remains a mys- 
tery '5| . These events have been seen by the Akeno Giant 
Air Shower Array (AGASA) 3J, Fly's Eye [Ij, Haverah 
Park HiRes and Yakutsk experiments. In- 
teraction of charged as well as neutral primaries with the 
cosmic background photons through photo-production of 
pious, photo-pair production, and inverse Compton scat- 
tering leads to a severe energy degradation and con- 
strains the source of a ~ 10^° eV charged particle to dis- 
tances less than ~ 50 Mpc (the high energy GZK zone). 
The lack of the most commonly hypothesized accelera- 
tors viz. accretion disks and jets associated with active 
galactic nuclei within the GZK zone has led to numer- 
ous alternate proposals for the origin of the Ultra-High 
Energy Cosmic Rays (UHECRs). 

These explanations fall into two broad categories of 
top-down and bottom- up scenarios. Top-down models 
generally rely on the decays of very heavy particles, usu- 
ally remnants from the early universe |3- Bottom- up 
models are based on the acceleration of normal charged 
particles to extremely high energies in astrophysical ob- 
jects not otherwise known to possess such effective accel- 
eration mechanisms @. These accelerators are expected 
to show a correlation with some known class of astro- 
physical objects in the local universe. 
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Supernovae and their remnants have long been sus- 
pected as the source of ordinary cosmic rays below 
4 X 10^^ eV. Two models based on gamma ray bursts 
(GRBs) and magnetars for the acceleration of UHECRs 
have been proposed that attribute these particles to brief 
bursts of particle emission associated with the forma- 
tion of compact objects in core collapse supernovae. In 
the GRB model, Milgrom and Usov Vietri |T3 and 
Waxman independently suggested that the relativis- 
tic shock waves associated with GRBs might be the ac- 
celeration sites for these cosmic rays. The relativistic 
shock model has been useful in accounting for the phe- 
nomenology of GRBs . Recent observations of GRB 
afterglows have supplied strong evidence that at least the 
long burst sources are associated with supernovae occur- 
ring in external galaxies 'l4'| . The prompt association of 
the afterglow of GRB030329 with SN2003dh gives strong 
support to the prompt formation of the collapsed object 
that drives the outburst, either a black hole, as in the 
coUapsar model IlJ j or a rapidly rotating magnetar ; 
the latter model works only if the jet that drives the 
GRB contains only part of the rotational energy deliv- 
ered, with the rest driving the supernova envelope, as in 
the magnetar model for UHECRs proposed in Ref. [l7| . 

Arons 0| showed that magnetars, a sub-class of neu- 
tron stars with ultra strong surface magnetic fields [B ~ 
10^^ Gauss) born in core collapse supernovae and hy- 
pothesized to have initial spin rates close to the centrifu- 
gal breakup limit, are unipolar inductors that can accel- 
erate particles in their relativistic winds up to energies 
~ 10^^ eV. They have an electrodynamically defined in- 
jection rate sufficient to yield a quantitatively acceptable 
account of the observed UHE spectrum, if such events 
occur in external galaxies with core collapse supernovae. 
Blasi, Epstein, and Olinto 0| presented a galactic pulsar 
variant of this idea, and Milgrom and Usov ^3 earlier 
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showed that such objects have voltages (including volt- 
age drops in the wind) large enough to account for the 
maximum energies seen in UHECR. Since some of the 
magnetars in our gal axy are physically associated with 
supernova remnants [l9j, the association of magnetars 
with core collapse supernovae is plausible. 

Since both the GRB and magnetar models attribute 
UHECRs to compact objects formed in core collapse su- 
pernovae (probably in the Ib/c sub-class), the rate of 
events that give rise to UHECRs in a galaxy should be 
proportional to the supernova rate in galaxies that have 
core collapse supernovae. Likewise, the core collapse su- 
pernova rate has been observed to be strongly correlated 
to the star formation rate (SFR) [23|, a correlation read- 
ily understood from the fact that such supernovae have 
short lived, massive stellar progenitors. This correlation 
leads us to explore the connection between the UHECR 
events and infrared bright spirals that are known to be 
active star forming regions. 

In this paper we test the hypothesis that the origin 
of the UHECR events lies in reasonably prompt accel- 
eration following a core collapse supernova, by study- 
ing the correlations between the angular locations of 
UHECR events on the sky and the positions of galaxies 
with a high SFR using the two-dimensional Kolmogorov- 
Smirnov (KS) test. We choose the IRAS Point Source 
Catalogue (PSCz) for our galaxy sample 0. The cata- 
logue consists of 14800 galaxies with measured redshifts 
and far infrared luminosities above a flux limit of 0.6 
Jansky at 60/im. The sky coverage is ~ 84%. 

Seciniprovides a summary of the acceleration theories. 
Sec mil discusses observational and theoretical evidence 
for the expected formation rate of the relevant compact 
objects in infrared bright spirals, and discusses the use 
the IRAS PSCz catalogue and the AGASA sample in 
our study. Sec IV describes the statistical methods and 
results on the data and theoretical models. Sec V and VI 
provide further discussion and conclusion of our results. 



II. COMPACT OBJECT PARTICLE 
ACCELERATION THEORY 

A. Gamma Ray Bursts 

Gamma ray bursts, and especially their afterglow emis- 
sion, clearly reflect the excitation of the magnetized 
gaseous medium surrounding an explosive release of en- 
ergy from a compact object of at least stellar mass [T^j . 
The observed relativistic expansion rates {e.g., |22j |') sup- 
port that release having given rise to relativistic outflow, 
well modeled in the afterglow phase as a relativistic blast 
wave propagating through a surrounding circumstellar or 
interstellar medium. In analogy to the apparent ability of 
ordinary supernova remnant shocks to accelerate a small 
fraction of the particles in the surrounding medium to 
relativistic energies, the authors in Refs. E3j E3 

suppose that relativistic shocks, either that of an external 



blast wave expanding with Lorentz factor 7 ~ 100— 1000 
which creates the afterglow or those of the internal shocks 
with 7 ^ 2 — 10 associated with the prompt GRB, are 
responsible for the acceleration of the UHECR. The lumi- 
nosity density for UHECR is [/uhe 5 x lO"^"* ergs/Mpc^- 
sec 0, |23|, about a factor of 5 greater than the pho- 
ton luminosity density of GRBs (~ 10*"* ergs/Mpc'^-sec 
23] but well below the mechanical luminosity density of 
GRBs, ^ 10'*^ ergs/Mpc'^-sec, assuming ~ 1% efficiency 
in converting internal shock energy into photon energy 
[2^ ]. The synchrotron model of GRB emission provides 
the only evidence for particle energies in a GRB or its 
afterglow large enough to correspond to ions being at 
relativistic energies. 

The energetics of this model are demanding but not ex- 
orbitant. The total observed GRB rate is n^^^i'GRB/b ~ 
0.5 X 10~^ Mpc"'^ yr~^ 25], where n^^^ is the num- 
ber density of galaxies hosting GRBs, i^grb is the GRB 
rate per host galaxy, and /b = f2b/47r is the frac- 
tion of the outflow beamed into solid angle f2b, with 
/b ^ 0.01 inferred in many afterglows. Lon g bu rsts pro- 
vide about 70% of the observed GRB rate |20], that is, 
nl^VcRB 0.35 X 10^9 Mpc-3 yr-i/Jj,, where n°«^ 
is the number density of galaxies with core collapse su- 
pernovae n^^° , and i^grb is some fraction of the overall 
core collapse supernova rate per galaxy. Then the UHE 
energy demanded of each burst source, if long bursts are 
the sole contributors, is AiJ^yf = f/uHE/"g'^''^'GRB = 
10^^/b ergs = 10^^(/b/0.01) ergs. This energy is some- 
what in excess of the energy put into photons in each 
burst. As is obvious from the energy budget of rotation 
powered pulsars and, to a lesser extent, of AGN with well 
developed jets, there is no fundamental difficulty with 
particle acceleration power exceeding photon luminosity 
in a compact object, especially if the outburst draws its 
power from the pressure of large scale electromagnetic 
fields. 

The microphysics arguments advanced for GRBs be- 
ing a source of UHECR depends on the theory of parti- 
cle acceleration in shocks. Most authors have employed 
the much studied theory of diffusive Fermi acceleration 
in shocks in the test particle limit to estimate maxi- 
mum energies and spectra, with the efficiency of shock 
acceleration left as a free parameter. Here, we reiterate 
the well known remark p.7| that the maximum energy a 
charged particle can attain in a relativistic shock wave is 
~ ZeBRs — 3 X 10^*'i3ioi?si7 eV, where Z is the atomic 
number of the charged particle, e is the electronic charge, 
-Bio is the systematic magnetic field of the medium into 
which the shock propagates in units of 10 Gauss, and 
i?si7 is the radius of the shock wave in units of 10^^ 
cm; i?si7 ^ 1 is typical of the shock size during after- 
glow emission. This energy corresponds to a particle's 
Larmor radius having reached the size of the shock, at 
which point the particle can move through the motional 
electric field as a freely accelerating particle and expe- 
rience the full electric potential drop. Acceleration at 
internal shocks, which occur at much smaller radii, re- 
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quire concomitantly larger magnetic fields, if GRBs are 
to be the acceleration sites for the highest energy parti- 
cles. Such fields are much larger than found anywhere 
in the interstellar media of galaxies, forcing the propo- 
nents of this model to appeal to rclativistic shock prop- 
agation in unusual environments, such as the B oc r^^ 
field in the possibly strongly magnetized circumstcUar 
winds which might surround a presupernova star 28] or 
the similarly structured field in the relativistic wind from 
a newly formed neutron star (Ref. [23| and references 
therein) or from a transient rnagnetized disk around a 
newly formed black hole (Ref. [13 )• 

Granted the existence of such environments, the par- 
ticle spectrum created by such a shock may be a power 
law with f{E) = E~p, p = 2 — 2.5; detailed studies of 
the test particle problem yield p — 2.2 — 2.3 2S]- Mod- 
els of the acceleration physics which incorporate some 
information about the magnetic turbulence that scatters 
downstream particles back across the shock front indi- 
cate possibly steeper spectra |23- Such turbulence is an 
essential ingredient if the mechanism as applied to shocks 
with flow Lorentz factors less than 10'^ are to accelerate 
particles to UHECR energies. The efficiency with which 
such shocks can convert flow energy into accelerated par- 
ticle spectra is unknown, although judging from the ac- 
celeration efficiencies for electrons and positrons found in 
the much more highly relativistic shocks terminating the 
winds from rotation powered pulsars, efficiencies as high 
as 1-10% appear possible [3lj . 



B. Magnetars 

Magnetars are a special class of neutron stars which 
form in core collapse supernovae. Studies of their associ- 
ation with supernova remnants in our own galaxy suggest 
the overall birth rate of magentars is on the order of 1- 
10% of the total core collapse supernova rate To"]; this 
rate is comparable to the rate of type Ib/c supernovae. 
Therefore, the magnetar birth rate is in effect propor- 
tional to the supernova rate in a galaxy. The theory 
described in 0| suggests injection of charged particles 
with maximum energy 



Err 



Ze$i = Ze 



3 X lO^^ZBi^n'l eV. (1) 



Here is a magnetar's surface magnetic field, Bi^ — 
B^/IQ^^ Gauss, fij — 10'*ri4 is the initial angular ve- 
locity of the neutron star, i?* is the stellar radius, and 
c is the speed of light. The initial rotation period is 
Pi = 0.64/ri4 msec. The ions actually gain their energy 
in the relativistic wind electromagnetically expelled from 
the neutron star at distances, r, much larger than the radii 
of the star and its magnetosphere, which allows them to 
avoid catastrophic radiation losses; the electric potential 
in the wind is rE = rB = If Z = 1 — 2, as may be 
suggested by the shower data 0| , one requires Pi < 2 — 3 



msec for the model to be viable. Studies of the termi- 
nation shocks of pulsar winds show that the relativistic 
ions gain their energy in the wind interior to the winds' 
termination shocks, rather than at that (reverse) shock, 
or at the blast wave driven into the surrounding inter- 
stellar medium 32]. The magnetar model in ^3] adopts 
this acceleration site by analogy to pulsars, rather than 
appealing to the relativistic blast wave driven by a newly 
born magnetar's fields into the surrounding interstellar 
medium. 

The star spins down rapidly, under the influence of 
electromagnetic and gravitational wave torques. As it 
spins down, the voltage and the particle energies de- 
cline. Summing over the formation and spindown event, 
one finds a per event injection spectrum proportional to 
f{E) = E-^[1 + {E/Eg)]~^ for E < E^^^. Here E^ mea- 
sures the importance of gravitational wave losses (calcu- 
lated for a star with static non-axisymmetric quadrupole 
asymmetry) in spinning the star down. These losses have 
significance only at the largest rotation rates, if they 
matter at all. When they exert torques larger than the 
electromagnetic torque, the star spends less time at the 
fastest rotation rates and therefore spends less time ac- 
celerating the highest energy particles, thus causing a 
steepening in the spectral slope at the highest energies. 
If the star has an internal magnetic field even stronger 
than the already large surface field, equatorial elliptic- 
ities ec in excess of 10^"^ can exist, in which case Eg 
would be less than E'max- We consider three cases of 
energy loss due to gravitational radiation (GR) in the 
model: no GR loss (eo = 0, i?g = oo); moderate GR 
loss (ec = 0.01, = 3 X 10^" eV); strong GR loss 
(ee = 0.1, = 3 X 10^** eV). 

The total number of particles injected per event is 



ZeB^, 



10^ 
ZBi5 
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where we have assumed a stellar radius of 10 km and 
a moment of inertia / — 10^^ cgs. In this respect, the 
magnetar model is much more specific than the GRB 
model, since the ion flux is uniquely flxed by the electro- 
dynamics of the unipolar inductor. In order to tap the 
motional EMF of the wind, the particles must cross mag- 
netic field lines and move parallel to the electric field. In 
Ref. 0, it was suggested that the current sheets in the 
structured wind are unstable to the formation of large 
amplitude electromagnetic waves, whose ponderomotive 
force drives the particles across -Bwind, thus causing the 
energy gain, a mechanism analogous to those studied in 
the design of advanced accelerators and of some interest 
for ion acceleration in the observationally much better 
constrained pulsar winds. 

The rate at which galaxies inject UHECR into the uni- 
verse in this model then is h„ — i^migM^gaiaxy, where 
"galaxy « 0.02 Mpc"^ [sj , Ni is givcu by Eq. ©, and 
'^mag the birth rate of rapidly rotating magnetars per 
galaxy. If intergalactic propagation is unaffected by scat- 
tering in an intergalactic magnetic field, multiplication 
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of the source spectrum q{E) oc ncrf{E) by the energy 
dependent GZK loss time yields a spectrum received at 
the Earth in reasonable accord with the existing obser- 
vations of UHECR, as shown in Figure 1 of Ref. [l7l| . 
if t'^lg ~ 10^^ yr~^- That fast magnetar birth rate lies 
between 1% and 10% of the total magnetar birth rate 
inferred for our galaxy, and about 0.1% of the total core 
collapse supernova rate in average star forming galaxy, 
- 10-2yr-i H. 

In the fast magnetar event rate was applied to 
all normal galaxies, whose space density is ng = 0.02 
Mpc~^ [ill, yielding an event rate ngi/^"* ~ 2 x 10"'^ 
Mpc~^-yr~^. The luminosity density of injected energy 
from the wind, which is about equal to the energy ex- 
pended in accelerating the ions, is 10^^, 0.9 x 10*^, and 
0.9 X 10**^ ergs/Mpc^-yr^^, in the no GR, moderate GR 
and strong GR cases, assuming the same i'^^^ in all three 
applications of the model. In fact, since v^^^^ is an other- 
wise unknown fraction of the poorly known i/m, the value 
of v^^^ was adjusted in each case so as to provide the 
best agreement between the model's UHECR spectrum 
and the observations. The UHECR energy per event 
is much smaller than is required for the GRB model, 
whose effective rate of occurrence (including beaming) 
is much less than for magnetars. However, there are 
many normal galaxies within every resolution element 
of the arrays used to accumulate the UHECR events, 
making angular correlations between all normal galaxies 
and UHECR events difficult in testing the magnetar (and 
the long burst GRB) hypothesis for the UHECR's origin. 
The luminous infrared galaxies offer a better prospect for 
testing these models, since they are rarer than ordinary 
galaxies but still contribute a substantial fraction of the 
metagalactic supernova rate (see below). 



III. STUDY SAMPLE AND SKY COVERAGE 

In this section we discuss how the far infrared luminos- 
ity of a galaxy serves as an indicator of its core collapse 
supernova rate and why the IRAS catalogue is chosen as 
the galaxy sample. We then discuss the use of AGASA 
events as the UHECR sample. 



A. Infrared Galaxies 

Type Ib/c and type II supernovae are considered good 
indicators of the SFR in a galaxy. The progenitors of 
these core collapse supernovae are young, massive stars 
that are abundant in active star forming regions. Ob- 
servational support for relating the SFR to the far in- 
frared luminosity (A > 20/im) has come from study of 
the IRAS catalogue of ^ 15000 infrared bright galaxies. 
The method relies on the fact that a large percentage of 
the bolometric luminosity of young stars is absorbed by 
interstellar dust and then emitted as thermal radiation 
from the warmed dust in the far infrared spectrum [20| . 



The absorption cross section of dust is strongly peaked in 
the ultraviolet, where the young stars dominate the ra- 
diation. This makes the thermal far infrared luminosity 
of a galaxy a good proxy for the SFR. 

Specifically, we take the supernova rate and therefore 
the birth rate of stellar mass collapsed objects to be pro- 
portional to the luminosity of a galaxy at 60/im. This 
wavelength is chosen to reduce the contamination from 
the radiation field of older stars at ^ lOfim as well as the 
cirrus emission at > lOOyitm due to stellar optical/near IR 
radiative excitation of large molecules in the interstellar 
medium |20l |. 

Luminous infrared galaxies with high SFRs and super- 
nova rates offer the best chance of finding a correlation 
between UHECR event source directions and sites of core 
collapse supernovae. We use the IRAS PSCz catalogue 
01 as the host galaxies for GRBs and magnetars. This 
catalogue has a nearly uniform coverage of ~ 84% of the 
sky, where most of the unmapped area lies within ~ 10° 
north and south of the galactic plane. The sample pro- 
vides redshifts for all galaxies above a flux limit of 0.6 
Jansky at 60/im. 

These galaxies contribute a substantial fraction of the 
overall core collapse supernova rate. From the IRAS cat- 
alogue, these objects have space density dnf^j-g/dLio w 
10-3L-i±°-^ Mpc-3, where Lm = Lfir/IO^^ Lq and 
Lfir is the far infrared luminosity |35[. The integral 
number of these galaxies (1 < Lio < 10'^) is rifirg w 0.005 
Mpc~^, which should be compared to the total space den- 
sity of galaxies, rig « 0.02 Mpc~^ [S^; most of these are 
spirals and other galaxies with some degree of ongoing 
star formation. 

The general core collapse supernova rate per galaxy, 
summed over all galaxy types, is ~ 0.011 SN per 
galaxy- year j34]. which yields a volume averaged rate of 
supernovae of Ssnng « 2.2 x 10"'' SN/Mpc^-year. The 
supernova rate for galaxies with far infrared emission is 
Sf^l = 2.5 X lO^Uio SN/(FIR galaxy)-year |3|. Inte- 
grating this supernova rate over the FIR galaxy luminos- 
ity function yields S'f^firg « 0.7 x lO-^K^hsciLTo'''' /300) 
SN/Mpc^-year. -ftTobsc, the correction for supernovae 
missed in the existing optical and near infrared supernova 
detection surveys, might be as large as 10, and probably 
is at least as large as 3 |0. Thus, the luminous infrared 
galaxies contribute at least 28% (i^obsc = 1) of the total 
supernova rate, with a total space density only 25% that 
of all normal galaxies. Furthermore, since the brightest 
infrared galaxies (-Lfir > 10^^ Lq) dominate the contri- 
bution from all FIR galaxies, and these are quite rare, 
with space density ~ 4 x 10~* Mpc~'^ [s^, correlations 
of UHECR arrival directions with the sky positions of 
the luminous IRAS galaxies offers a promising opportu- 
nity to test the hypothesis that UHECR acceleration has 
something to do with core collapse supernovae, as is im- 
plied by the GRB shock and magnetar unipolar inductor 
models for the acceleration sites. 
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B. UHECR Events 

For the UHECR events, we choose the AGASA experi- 
ment because among the ground based arrays, this exper- 
iment has the best angular resolution (1.8° FWHM js^) 
and provides the largest single dataset above 4 x lO""^^ eV. 
The current sample of AGASA UHECR events contains 
60 events above 4 x 10^^ eV, among which 57 events are 
published in Ref. 38] and 3 recent events with energies 
greater than 1 x 10^° eV are listed at the AGASA website 

In order to test whether the AGASA events are con- 
sistent with their origin being from IRAS galaxies, we 
must first take into account the different sky coverage by 
AGASA vs. IRAS and the possible variations in detec- 
tor exposure efficiency in AGASA. Continuously operat- 
ing ground-based arrays such as AGASA typically have 
a uniform coverage in sidereal time and therefore little 
exposure variations in right ascension (RA) |37j . The de- 
tector exposure in declination, however, depends on the 
latitude of the experiment and drops to zero beyond 
certain maximum angle, ^max- For AGASA, = 35.8° 
and ^max = 45°. Beyond 45°, the energy determina- 
tion of the events becomes uncertain. In terms of these 
parameters, the detector efficiency as a function of dec- 
lination, 9, is given by j40ll4]|. 



de 



oc (cos 6*0 cos sin + am sin 9osin9) , (3) 



where am is given by 




and 



if ^ > 1 
if ^ < -1 
otherwise 



^ _ cos flmax — sin 6*0 sin ( 
cos 9o cos 9 



(4) 



(5) 



For comparison, we find the observed declination dis- 
tribution published in Ref. [s^l for the number of AGASA 
events above 1 x 10"^^ eV to be well fit (to ~ 1%) by 

dNo _ ( -{9 + 10){9 - 80)/35.5 , if - 10° < 6i < 80° 
d9 ^ [ , otherwise 

(6) 

We have verified that the declination distribution for 
events above 4 x 10^^ eV is similar to Eq. ©. This 
is consistent with the triggering efficiency being uniform 
above 1 x 10^^ eV and the UHECR events not showing 
any appreciable large scale anisotropy 37, 42, 43|. We 
use the theoretical window function in Eq. ||2J) for the 
statistical tests in Sec. IV and discuss the implications of 
using the observed declination distribution in Sec. V. 

Fig. 1 shows the sky distribution of the AGASA 
UHECR events and the area not mapped in the PSCz 



catalogue. Six AGASA events lie in these unmapped re- 
gions and will therefore be excluded, leaving 54 AGASA 
data points above 4 x 10^^ eV for our study. 



IV. STATISTICAL TESTS 

We use the two-dimensional KS test [i^. l45l | to test 
whether the AGASA UHECR events above 4 x lO^^ eV 
are consistent with the cosmic-ray intensity distribution 
predicted by the GRB and magnetar models in infrared 
bright spirals discussed in Sees. II and III. In this sec- 
tion we first describe our calculations of the expected 
UHECR distribution from the models and then discuss 
the statistical test and the results. 



A. Model Predictions 

The spectrum of cosmic rays at Earth depends on both 
the injection spectrum at the acceleration sites and pos- 
sible particle interactions that would degrade the ener- 
gies of the primaries between the source and Earth. As 
discussed in Sec. II, we model the injection spectrum to 
be f{Ei) = ErP, p = 2-2.5 for the GRB models and 
f{Ei) = Er'^[l + {Ei/Eg)]-^ for the magnetar models. 
For the latter we will explore the three gravitational radi- 
ation cases described in Ref. [l3|: no GR loss {Eg oo), 
moderate GR loss {Eg — 3x 10^° eV), and strong GR loss 
{Eg = 3 X 10^* eV). To estimate the energy degradation 
during propagation, we use the probability Pp{r,Ei]E) 
that a proton created at distance r with energy Ei would 
arrive at Earth with an energy greater than E. It mea- 
sures the proton energy degradation due to the GZK loss 
processes discussed in Sec. I. The same function was used 
in our earlier paper (461] and has been calculated using 
Monte Carlo simulations ^3] for the parameters consid- 
ered here. 

We combine the injection spectrum f{Ei) and propa- 
gation loss Pp, and calculate the cosmic ray flux at Earth 
[in (eV m^ s sr)"-"-] by 



F{E,n) 



K{^) / dEi / dr[l + z(r)] = 
^0 Jo 
dPp{r,Ei;E) 



xf{EO 



dE 



(7) 



Here z is redshift, -Emax is the maximum proton injection 
energy, and K is a normalization constant which depends 
on the cosmic ray injection model and the viewing direc- 
tion in the sky, fl. The KS test (see next subsection) 
depends only on the relative spatial distribution of the 
UHECR intensity in the sky, we therefore need not spec- 
ify K in the theories. 

Fig. 2 compares the cosmic ray spectrum F{E)E^ for 
various GRB and magnetar models. For illustrative pur- 
poses, the figure is constructed assuming isotropic injec- 
tion with a uniform comoving density, i.e., K is indepen- 
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dent of direction on the sky and is proportional to the 
average source density and the event rate per unit time 
at the present epoch. If UHECRs come from discrete 
sources such as galaxies, this procedure assigns the same 
cosmic ray luminosity to each source, an assumption we 
will relax below. The curves in Fig. 2 have been nor- 
malized so that all model predictions have the same flux 
at 4 X IQi^ eV. The GZK feature at £; > 4 x lO^^ eV 
becomes more pronounced for steeper injection spectra. 
(We note that Fig. 2 corrects a plotting error in the low 
energy spectrum {E < 10^^ eV) of the strong GR case in 
Fig. 1 of Ref. ItI.) 

Fig. 3 illustrates the expected dependence of the aver- 
age UHECR flux /J'""'' F{E)dE on the luminosity dis- 
tance, dL = cH^^ {1 + z)[z- z'^{l + qii)/2 + 0{z'^)], to the 
host objects for various models: GRBs with a E^"^-^ in- 
jection spectrum (long dashed), magnetars with no GR 
loss (solid and dotted), and magnetars with moderate 
GR loss (short dashed). Two lower energy limits -Emin 
are shown for comparison for the no GR loss model. For 
E^i-a = 10^" eV (dotted), the integrated flux falls off 
sharply above d ~ 50 Mpc due to the GZK loss pro- 
cesses. For i?min = 4 X 10^^ cV (solid), the cutoff is more 
gradual and occurs at a larger distance {d ^ 200 Mpc). 

In order to compare the distribution of UHECR events 
with the hypothesized discrete sources, we now need to 
provide an estimate of each IRAS galaxy's contribution 
to the UHECR flux. We do this by reinterpreting Eq. O 
to be the flux expected from each galaxy in the IRAS cat- 
alog, with the galaxies chosen to lie within the maximum 
GZK zone indicated by Fig. 3. That is, K in Eq. Q 
is a number chosen to be proportional to the 60 /xm lu- 
minosity for each galaxy with a catalogued sky position, 
all within the distance i?max, while the integral over r in 
Eq. Q is carried out for each galaxy separately, from 
to the luminosity distance di^{z) of each galaxy, using the 
catalogued redshift, z, for that galaxy. This is done for 
each of the three model injection spectra of the magnetar 
model, and for two possible GRB injection spectra shown 
in Fig. 2. Finally, we smooth the galaxy model's predic- 
tions of the UHECR flux as a function of sky position to 
the angular resolution of the AG AS A experiment, using 
a top hat angular smoothing profile. 

We include all IRAS PSCz galaxies within i?inax = 1 
Gpc, although our results are not sensitive to this up- 
per cutoff as long as it is greater than the GZK zone 
of 200 Mpc. We estimate Lfir by using the 60 ^m 
flux and the flux-luminosity relationship L — Aird^fpm. 
We use matter density rJm = 0.35, cosmological constant 
= 0.65, and Hq = 70kms"^Mpc"^ The luminosity 
distance, c^l differs from linear d = cz/Hq by ^ 15% at 
1000 Mpc and by ~ 3% at 200 Mpc. Our final results 
depend very weakly (~ 1%) on the exact cosmological 
parameters assumed in the calculations. 

Fig. 4a shows a linearly spaced contour map of the 
UHECR integrated flux predicted by our infrared galaxy 
source model computed from IRAS galaxies for a magne- 
tar model (no GR case). The SFRs (and hence magnetar 



rates) associated with the galaxies are assumed to be lin- 
ear in L-piR. For comparison, all 60 AGASA UHECR 
events above 4 x 10^^ eV are superposed. Fig. 4b shows 
the same quantity for a GRB model with a E~^ injec- 
tion spectrum. No corrections for the limited sky cov- 
erage and non-uniform detector efRciency in declination 
(see Sec HI) have yet been applied. 

The question at hand is, do the model and observed 
angular distributions of UHECR event rate on the sky 
significantly differ? That is, can we disprove the null hy- 
pothesis, that the AGASA events are drawn from sources 
located in IRAS galaxies? 

B. Kolmogorov-Smirnov (KS) Test 

Since the AGASA events do not obviously point to- 
wards individual galaxies, a statistical test is needed to 
quantify if the AGASA events are inconsistent with the 
theoretical UHECR distributions shown in maps such as 
Fig. 4. 

The two-dimensional KS test is a generalization of the 
much studied one-dimensional KS test. The latter has 
been discussed extensively in the literature and analytical 
proofs for its validity can be found in R is based 

on comparing the cumulative probability distributions of 
the sample H{x) and the model G{x). If Z?id is the 
maximum absolute value of H{x) — G{x) for all x, then 
the probability P that H{x) and G{x) are drawn from 
the same underlying distribution is given by 

lim P(> V^Du) = QiV^Did) (8) 

n— >oo 

and 



Q(A) = 2^(-l)^-le-2^"^^ (9) 

k=l 

which is valid if the number of data points in the sam- 
ple, n, is > 80. Monte Carlo methods must be used to 
tabulate the probabilities for smaller n. 

Note that a small computed value of P indicates that 
the null hypothesis is false, that is, the data did not come 
from a distribution represented by the model. Generally, 
P < 0.1 is considered reasonable proof that the null hy- 
pothesis is false. Larger values of P simply provide ev- 
idence that the data and model are consistent; the null 
hypothesis has not been proved false. 

The generalization of the KS test to two dimensions is 
not obvious because there is no unique cumulative prob- 
ability distribution in more than one dimension. At least 
two different generalizations of the test to two dimensions 
have been proposed 0, . Both methods have been 
verified with Monte-Carlo simulations although no ana- 
lytical proof for their validity exists. In this paper we use 
the probabilities tabulated in Ref. [l^. The test com- 
putes the difference between the fraction of data points 
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and the fraction of theoretical points in each of the four 
quadrants of the plane of the sky centered at a data point 
of angular coordinates {xi,yi), where the quadrants are 
defined by {x > Xi,y > yi), {x < Xi,y > yi), (x < Xi,y < 
yi),{x > Xi,y < yi). This procedure is repeated for all 
data points, i = 1, n, resulting in 4 x rt numbers for the 
differences between the data and model. The maximum 
absolute difference value is defined as I?2d- Monte-Carlo 
simulations are used to obtain the corresponding proba- 
bility P(> ^JnD-id) for a given and n 01- It is also 
shown in Ref. [4JI that the results are almost indepen- 
dent of the underlying distribution function. For large n, 
the probability has an analytic form i49l] : 

P(> ^D,,) = Q { ^ ^^'^ \lO) 

where Q(\) is given by Eq. ^ and r^c is the linear 
correlation coefficient of the 2-dimensional data points 
ixi,yi),i = l,n: 

r = - x)iyi - y) 

v/Ei(^i-s)VEi(2/i-y)'" 

where x and y are arithmetic means, and — 1 < < 1- 
The value rcc = 1 is for perfectly correlated x and y, 
which as expected, reduces Eq. (|1U|) to the same form 
as the one-dimensional expression in Eq. The value 
Tec = is for no correlation between x and y, as is the 
case for the AGASA events. 

In this study, we take the data points to be the an- 
gular coordinates {x,y) =(RA,dec) for the 54 AGASA 
UHECR events above 4 x 10^^ eV that lie within the sky 
coverage of the IRAS sample. The model predictions are 
the UHECR intensity calculated from IRAS galaxies as 
discussed in Sec IVA, and further weighted by the func- 
tion dNt/dO in Eq. Q to account for the non- uniform 
exposure of the AGASA detector in declination 6. We 
determine the KS parameter D2d by searching through 
the four quadrants centered at the coordinates (RA,dec) 
of each AGASA event. For each quadrant, we compute 
the difference between the fraction of the UHECR events 
that lie in that quadrant, and the fraction of the UHECR 
intensity predicted by the models from IRAS galaxies in 
the same quadrant. We then assign the maximal differ- 
ence from all 4 x 54 quadrants to be I?2d- We determine 
the probability P(> y/nD2d) from D2d by interpolating 
Table I of Ref. |4J| . We do not use the analytic expres- 
sions in Eqs . © and (|10|l that are adopted in the algo- 
rithm of |43 because they are accurate only for n ^ 20 
and P ^ 0.20. (For instance, we find a probability of 
32.9% using the subroutines provided in ^3 vs. 50.6 % 
from 44] for the magnetar-no GR model discussed be- 
low.) 

Our results from the KS test for various magnetar and 
GRB injection spectra for UHECRs are summarized in 
Table I. The SFR of each IRAS galaxy is assumed to 
be proportional to its Lfir- The consistency between 



models and data above 4 x 10 eV is generally high with 
KS probabihties in the range of 50.6% to 51.7% (third 
column in Table I); that is, the null hypothesis, that the 
angular coordinates of the AGASA events come from the 
set of angular coordinates of the IRAS galaxies within the 
GZK zone, is not disproved. If the data set is restricted 
to the 8 highest energy events above 10^° eV (out of 
11 total; 3 lie outside the IRAS coverage), however, the 
KS probabilities dive to between 0.30% and 0.45% (fifth 
column); that is, the models are rejected at > 99.5% 
significance level. 

The assumption of SFR cx Lfir is reasonable but is 
subject to some theoretical and observational uncertain- 
ties (113 and 0). We therefore test if the KS probabili- 
ties are sensitive to the linear relation of SFR to Xfir by 
examining some extreme cases where SFR cx ipiRi with 
a — 0,0.5,1,1.5, and 2. Table II summarizes the re- 
sults for a magnetar (no GR loss) and GRB {E~'^) injec- 
tion spectrum. The test shows only small variations for 
a ~ Q, 0.5, and 1, with the probabilities dropping sharply 
only for very steep power laws to 24.1% for a = 1.5 and 
0.56% for a = 2 for the magnetar (no GR loss) injec- 
tion spectrum. Similar results are found for the GRB 
injection spectrum. We conclude from this test that our 
results are robust to all reasonable range of a 1. 

Starburst galaxies are generally associated with lumi- 
nous far-infrared galaxies (Lfir ^ lO^^i©) [HI. To 
explore a connection between these galaxies and the 
UHECR events, we separate all IRAS galaxies within 
200 Mpc into two samples with LpiR above and below 
10"Lo. We find the KS probability to be 44.3% and 
60.1% for the above and below IO^^Lq samples, respec- 
tively. Starburst galaxies are therefore not ruled out as 
possible UHECR sites, although the KS probability does 
not favor starbursts more than ordinary IRAS galax- 
ies. Combining the two samples gives a probability of 
49.6%, which is very close to 50.6% discussed above for 
the galaxy sample out to 1 Gpc. This is to be expected 
as most contribution to the UHECR flux comes from the 
GZK zone as indicated in Fig. 3. 

We perform the KS test on a number of additional toy 
models vs. AGASA UHECR events to explore the dis- 
criminating power of the KS method on model assump- 
tions. In all cases we make sure to weight the theoreti- 
cal predictions with the same AGASA detector exposure 
function dNt/d9 for a fair comparison. In the first toy 
model, we assume an isotropic distribution of cosmic ray 
intensity and obtain a probabilities of 46.2% and 9.1% 
for the UHECR events above 4 x 10^^ and 10^° eV re- 
spectively to have been drawn from an isotropic distri- 
bution. In the second toy model, we ignore the redshift 
and LpiR of the IRAS galaxies and simply test their sky 
positions vs. the 54 AGASA events. We find a KS prob- 
ability of 54.5% for UHECR events above 4 x lO^^ eV 
and 10.9% for events above 10^° eV. In the third toy 
model, we take the sky positions of the 2702 GRBs in 
the Burst And Transient Source Experiment (BATSE) 
catalog (http://www.batse.msfc.nasa.gov/batse). Since 
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TABLE I: KS probabilities for UHECR events above 4 x 10^^ eV and above 10^" eV and the cosmic ray intensity from various 
models to be drawn from the same sample. 





E> 4x 10" eV 


£ > 4 X 10" eV 


E > 10^° eV 


E > 10^° eV 


Model 


-D2d(n = 54) 


P(> VnD2d) 


-D2d(n = 8) 


P(> VnD2d) 


Magnetar: no GR loss 


0.1482 


50.6 


0.626 


0.31 


Magnetar: moderate GR loss 


0.1472 


51.7 


0.615 


0.39 


Magnetar: strong GR loss 


0.1478 


51.0 


0.609 


0.45 


GRB: 


0.1475 


51.5 


0.626 


0.31 


GRB: E-'-^ 


0.1474 


51.6 


0.627 


0.30 



TABLE II: KS probabilities for UHECR events above 4 x 10^ eV and the expected cosmic ray intensity from different SFR 
models to be drawn from the same sample. 





Magnetar: no GR 


Magnetar: no GR 


GRB: E-^ 


GRB: E-'^ 


Model 


D2d(n = 54) 


P(> VnDid) 


D2i{n = 54) 


P{> V^D2d) 


SFRocL" 


0.145 


55.5 


0.143 


57.0 


SFRcxL''-'^ 


0.144 


55.7 


0.143 


57.0 


SFRcxL^° 


0.148 


50.6 


0.147 


51.5 


SFRocL^'5 


0.169 


24.1 


0.168 


26.3 


SFRcxL^" 


0.267 


0.56 


0.271 


0.48 



the sky coverage here is 100%, we use all 60 AGASA 
UHECR events in this test. We obtain a KS probabilities 
of 42.3%(E> 4 X IQi^ eV) and 5.2%(E> lO^" eV). A re- 
alistic calculation of the correlation between the BATSE 
GRBs and UHECR events from AGASA gives a negative 
result. (50j The numbers from these three toy models 
are consistent with the fact that IRAS galaxies, BATSE 
GRBs, and the AGASA events are all nearly isotropically 
distributed on the sky on large angular scales ^51^, .52J . 
We do emphasize that taking only the sky positions of 
source catalogs and correlating them with AGASA coor- 
dinates (as done in these toy models and occasionally in 
the literature) is naive since it ignores important energy 
loss processes of the primaries during their propagation 
and possible dependence of the UHECR flux on source 
properties such as their intrinsic luminosity. 



V. DISCUSSION 

In this paper we have found using the two-dimensional 
KS test that the AGASA data for the UHECR events 
with E > 4 X 10^^ eV are consistent (at ~ 50% sig- 
nificance level) with the expected cosmic ray intensity 
distribution in the sky predicted by models in which the 
acceleration mechanism for the primaries resides in GRBs 
or magnetars in infrared bright galaxies. When only the 
highest UHECR events ( E > lO^o eV) are used, how- 
ever, we find the KS probabihties to be < 0.5%, indicat- 
ing the models are rejected at > 99.5% significance level. 
When the energy cut on the UHECR sample is varied to 
£: > 5 x 1019 eV, > 6.3 x lO^^ eV, and > 8 x lO^^ eV, 



we find the corresponding KS probabilities to be 54.2%, 
6.73%, and 0.92% for the magnetar model (no GR loss), 
suggesting a better correlation of the UHECR events at 
lower energies. 

To establish whether the low KS probabilities at the 
highest energies are due to the paucity of data points or 
due to the failure of the models under consideration, we 
draw 100 random samples of 8 events each from the data 
in the energy range, 4x 10^^ eV < E < 10^" eV. There are 
46 events in this energy range in the region mapped by 
IRAS and the KS probability is 43.3% for the magnetar 
with no GR version of the model. For a fair comparison 
with the AGASA data, we consider only those samples 
for which the absolute value of the linear correlation co- 
efficient I Tec I defined in Eq. (|11() is less than 0.1. We find 
that the mean KS probability of the 100 random sam- 
ples is 52.8 % with a standard deviation of 24.4 %. It 
is therefore unlikely (at ~ 2a level) that the poor agree- 
ment between the data above 10^" eV and the models 
considered here is due to the small number of AGASA 
events above 10^'^ eV. The AUGER experiment .53] with 
a projected ~ 60 events per year above 10^*^ eV [5j| will 
help resolve whether these highest AGASA events indeed 
represent a 2(7 statistical fluke. Taken at face value, how- 
ever, the existing UHECR data suggest that at least the 
cosmic ray events above 10^" eV do not owe their origin 
to long burst GRBs, rapidly rotating magnetars, or any 
other events associated with core collapse supernovae. 

Our work can be contrasted with a previous work [H^ 
on the correlations between the UHECR events and lu- 
minous infrared galaxies, in which the authors used the 
two-dimensional generalization of the Smirnov-Cramer- 
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von Mises test. The conclusions of the two papers are 
in broad agreement, although our significance levels dif- 
fer from theirs partially because unlike the K-S test 
used here, the generalization of the Smirnov-Cramcr-von 
Mises test to two dimensions used in Ref. ,55j does not 
yield a unique significance level and there is a large scat- 
ter. Moreover, Ref. [s^ does not exclude the AGASA 
data points that lie in the unmapped region of the IRAS 
survey. To illustrate the importance of correcting for 
the detector exposure, we replace the window function 
in Eq. JS)) with a tophat function that is equal to 1 for 
-10° < 6* < 80° and zero otherwise, and find the KS 
probability to drop to 33.7% from 50.6% (for the magne- 
tar model with no GR loss). 

We also test the robustness of our results by using the 
observed window function given by Eq. (|SJ) instead of the 
theoretical function in Eq. We find this to have only 
a ~ 10% effect, changing the KS probabilities to 55.7%, 
57.5%, 45.2%, 42.8%, and 43.1% for the five injection 
spectra (from top down) listed in Table I. 

The significance of our comparison of the angular co- 
ordinates of the UHE events on the sky to the coordi- 
nates of the IRAS galaxies relies on rectilinear particle 
paths through intergalactic space. Could scattering in 
intergalactic magnetic fields affect our conclusion? An 
interesting constraint on the intergalactic field and the 
fluctuation correlation length from the GRB model is 
provided by the 8 doublet and 2 triplet UHECR events 
occurring with 3° of each other on the sky (3^ , among a 
total of 92 events with E > 4 x 10^^ eV that have been 
observed over the last 40 years. Since no such close pairs 
of GRBs have been observed, each UHECR doublet and 
triplet are therefore events from the single GRBs. How- 
ever, these events are separated by a decade or more, yet 
GRBs last dynamically only for a few months. The par- 
ticles from a GRB therefore must have undergone time 
delays > 10 years during their propagation through the 
intergalactic medium, an effect which requires some form 
of scattering along the propagation path. The same con- 
siderations apply to the magnetar model, although in 
contrast to GRBs, photon observations of magnetar for- 
mation do not exist, so we lack explicit observational 
evidence against multiple magnetars occurring within 10 
year time intervals at the same angular position. 

In addition, GRBs within the very high energy GZK 
zone are extremely rare. The GRB rate (long bursts only) 
is ^ 3.5 X 10"^*^ Mpc"'^ yr~^. Particles with energy above 
4x 10^^ eV must have been emitted within the GZK zone, 
D < Dgzk = lOOE^o^ Mpc (3 > £^20 > 0.4), assuming 
straight line propagation. Here E20 is the energy in units 
of 10^" eV. Then an average of ~ 700i?2o/10^ years pass 
between each separate GRB contribution to the flux of 
such extremely high energy particles observed at Earth. 
Therefore at the highest energies, the contributions from 
the very rare nearby bursts must be spread out in time, if 
the observed event rate is to be attributed to this source. 
The magnetar model, with its larger event rate, does not 
require temporal dispersion of each burst of particle in- 



jection in order to account for the events observed above 
10^*^ eV. Such dispersion is required, however, in order 
to account for the doublets and triplets, since each mag- 
netar injection event lasts only a few days 

The most likely origin for such time delays is 
charged particle motion in an intergalactic magnetic 
fleld of average strength, iJigm, and irregularities with 
amplitude SB Big^ and correlation length Ic- 

One can show, by simple random walk arguments, 
that since a particle is deflected by the small an- 
gle 59 w Ic/rh — (ZeBig^lc/E) ^ 1 in cross- 
ing a single region of correlated field, the parti- 
cles from an impulsive event at distance D have ar- 
rival angles smeared over an angular width O = 

(DUrl)'/^ = OA°l'j;iZBao/EUD/DozK)'/^. Here 
-Siio = Agin/10^^° Gauss and ^Mpc = ^c/1 Mpc. 
Simple geometry shows the time delay associated with 
this angular scattering to be ATd = {D/2c)Q^ = 
365{D/Dgzk)'^{Z'^B?^o/EIq)Impc years. At the same 
time, the scattering cannot smear the inferred angular 
position of the source by an amount substantially larger 
than the angular resolution of the AGASA array. The 
requirements 8 < 3° and ATd > 10 years constrain the 
intergalactic field and its fluctuation correlation length 
to lie in the range ^EI^{DgzkI Df/^ > ZBuqIU^c > 
OME^g{DGZK/D), 0.4 < E20 < 3, if an impulsive 
event model is to be a viable candidate for the origin 
of UHECR. See Ref. [56j for a complete analysis of this 
issue, including a discussion of likely values for the inter- 
galactic field and the effect of transport through galaxy 
clusters as well as the general intergalactic medium. 

To estimate how much magnetic fields can effect our 
results, we take all 60 AGASA data points above 4 x 10^^ 
eV and randomly displace them by 3°. After rejecting the 
6 events that fall in the unmapped region of IRAS, we 
find a KS probability of 46.3 % for the no GR magnetar 
model, a result consistent with our earlier finding. A sim- 
ilar analysis for the AGASA events above 10^° eV yields 
a KS probability of 0.58 % comparable to the earlier 0.31 
% obtained without taking into account any deflection 
due to magnetic flelds. These results suggest that our 
conclusions do not depend on the deflection of cosmic 
rays due to magnetic fields for the values discussed in 
the preceding paragraphs. 



VI. CONCLUSIONS 

The two dimensional KS test applied to the angular 
positions of the 60 AGASA events with energies above 
4 X 10^^ eV and the angular positions of the far infrared 
galaxies in the IRAS PCSz catalog show that the sources 
of these particles are consistent with impulsive events 
associated with core collapse supernovae. The KS prob- 
ability for this association is ~ 50%, too large to reject 
the hypothesis. Possible physical models of accelerators 
associated with such supernovae are relativistic shock ac- 
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celeration associated with long burst GRBs and electro- 
magnetic processes in the relativistic winds from newly 
born, rapidly rotating magnetars. 

However, application of the same test to the 11 
AGASA events with energy above 10^'^ eV yields a KS 
probability of < 0.5%, showing that the existing set of su- 
per GZK events are not consistent with the GRB model, 
the magnetar model, or any other model based on im- 
pulsive events associated with core collapse supernovae. 
At the ~ 2(7 level, this inconsistency is unlikely to have 
occurred by chance. However, 2tT results in high energy 
particle and photon astrophysics are notorious for their 
lack of robustness. Therefore, this highest energy result, 
which uses a very small number of cosmic ray events, is 
perhaps best regarded as a prediction that will be tested 
by the AUGER experiment. That experiment will yield 
enough events above 10^° eV to show whether the con- 



sistency between these compact object models and the 
UHE data using the whole set of results from 4 x 10^^ 
eV and above, or the inconsistency between models and 
data using only the super 10^" eV data, is the correct 
conclusion. 
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FIG. 1: Sky distribution of the 60 AGASA data points above 4 x 10^^ eV (filled circles) and the approximate areas not mapped 
by IRAS (black bands) . Six AGASA events lie in the unmapped IRAS region, which we exclude from our study. 
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FIG. 2: Predictions for the cosmic ray flux as a function of energy assuming an isotropic distribution of sources with uniform 
comoving density for the five different injection spectra analyzed in this work: magnetar models with no GR loss (solid), 
moderate GR loss (short dashed), and strong GR loss (dotted), and GRB models with E~'^ (long dashed), and E~^'^ (dot 
dashed). The data points (squares) are taken from Ref. for an earlier set of 30 AGASA UHECR events. The GZK 

suppression is more important for steeper injection spectra. 
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FIG. 3: The integrated flux as a function of the source distance for the no GR loss injection spectrum for E > 4 x 10 eV 
(sohd) and E > lO^'' eV (dotted). In the latter case the integrated flux falls off very sharply for d > 50 Mpc. Also shown are 
the integrated flux above i? > 4 x 10^^ eV for the moderate GR (short dashed) and E~'^'^ (long dashed) injection spectrum. 
The upper cutoff in the injection spectrum was taken to be 3 x 10^^ eV for all cases in accordance with [l7||. 




FIG. 4; Contour plot of the expected integrated cosrnic ray flux from IRAS galaxies smoothed to the AGASA angular resolution 
of 1.8°. The 60 AGASA UHECR events above 4 x 10^'' eV are shown as stars. Two models arc shown for comparison: the 
magnetar model with no GR loss (top) and the GRB model with injection spectrum. No corrections for the incomplete 
IRAS sky coverage and non-uniform AGASA detector efficiency in declination (see Sec III) have been applied. 



